**************Weather Effects on Forestry Production in Southern China*******************
****Table S4 ******************
	global raw "~/agriculture_temperature"
	use $raw/data/forest_weather_idw_pro_region, clear	


  *set temperature bin:<10℃，15-25℃，>25℃
  gen dd5bin_less_10=dd5binNeg+dd5bin0+dd5bin5
  gen dd5bin15_25=dd5bin15+dd5bin20
  gen dd5bin25_45=dd5bin25+dd5bin30+dd5bin35
  drop dd5binNeg dd5bin0 dd5bin5
  drop dd5bin10 dd5bin15 dd5bin20
  drop dd5bin25 dd5bin30 dd5bin35
  sum dd5bin*
	xtset pro year

  *South：dd5bin_less_10 dd5bin0_25 dd5bin25_45 
 	     estimates clear
		  reghdfe ln_timber_production dd5bin_less_10 dd5bin15_25 dd5bin25_45 ///
			c.pre##c.pre c.ssd##c.ssd c.rhu##c.rhu c.prs##c.prs c.win##c.win ///
			c.year##c.year ///
			if inrange(year,1995,2015) & region==2,  ///
			absorb(pro year) cluster(pro) 
			estimates store n_`x' 
esttab n_* using $raw/revision_result/south_products.rtf,replace  b(4) se(4) keep(dd5bin* *pre* *ssd* *rhu* *prs* *win* ) star(* 0.10 ** 0.05 *** 0.01) scalar(r2)
